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Abstract 

In recent years there has been some interest in applying ideas and methods taken from Physics 
in order to approach several challenging mathematical problems, particularly the Riemann Hy- 
pothesis, perhaps motived by the apparent inaccessibility to their solution from a full rigorous 
mathematical point of view. Most of these kind of contributions are suggested by some quantum 
statistical physics problems or by questions originated in chaos theory. In this note, starting from 
a very simple model of one-dimensional lattice gas and using the concept of equilibrium states as 
being described by Gibbs measures, we link classical statistical mechanics to the Riemann Hypoth- 
esis. 
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The link between physics and number theory [l(,[2j has been enriched in the few last 
years, particularly because of the possibility of enlighten, from a different point of view, 
some challenging mathematical problems, such as Riemann Hypothesis [3 ]• Most of the 
contributions in this sense come from the quantum side of statistical physics and chaos 
theory (see, for example, [4j and references therein). In this note we approach to the Riemann 
Hypothesis from a classical statistical mechanics model through the concept of equilibrium 
states as derived from a finite system variational principle. 

Physically our model is very simple: a one-dimensional lattice gas in the grand canonical 
ensemble. The lattice is an interval of the natural numbers: [1,M] C N where M is large 
enough (eventually, in the thermodynamic limit, M — > 00). The system is in contact 
with a particles reservoir characterized by the chemical potential fi and a heat reservoir at 
temperature T. For iV particles the configuration of the system is given by u; = (ojn, N) 
where uj^ = (h, h, 4 • • , «jv)- The coordinates i a (a = 1,2, ■ ■ ■ , N) take values in [1, M] and 
iV ranges, in principle, between and 00. The set of such configurations, the configuration 
space, will be denoted A. 

We assume that each site in the lattice can have at most one particle, so we are considering 
a hard point pair potential 

/. . x foo if i a = ig ... 
^0 if i a ^ ip. 

Also we assume that the particles are subject to a one-point potential such that they can 
just occupy sites in the lattice which are prime numbers: 

Ui u \ - J if *« is P rime (2) 
\oo if i a is composite. 

By the way, an example of explicit function with this form is|5|]: Ui(i a ) = 
— log[sin 2 |7T [(i a — l)!] 2 /i a } / sin 2 {7r/z a }]. The energy for the configuration u is u (u) = 

Mathematically, a state of the system is a probability vector v — {y (u) \u G A). The 
set of all states is denoted by M.. In the state v the system has the mean energy 
v ( u ) := SwgA v ( w ) u The (grand) partition function of u is defined S(/3,/x) : = 

^^ eA exp [— /3 (m (w) — yuA^)] where /3 = (fc^T) -1 with k B the Boltzmann constant. For 
the parameters /3, M and \x the Gibbs measure is defined 
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Given a measure v we have the entropy H (u) := — X] weA v ( w ) log 27 ( w ) an d the grand 
potential 

OH := v (u) - ifT (i/) - ju/ (iV) (4) 

withiz(iV) :=J2^ a u(uj)N. 

To introduce the notion of equilibrium state we consider a variational principle for finite 
systems (|A|=finite) according to which, for given energy u and parameters 0, M and /i, the 
Gibbs measure satisfies 

n[u ] = ~]ogE(^M lf i)= inf L{u)-\h{v)-^{N)\. (5) 

A measure that attains this infimum is called an equilibrium state. Gibbs measure vq is thus 
an equilibrium state. The result given by Eq.flSJ) is easily demonstrated [6] by using Jensen 
inequality applied to the concave function xH-lnx. The principle can be expressed saying 
that for any measure v is ^ Q [i/ ] = — jj log 5 (/?, M, /x) with equality if and only if 
v = u . 

Turning to our model, the grand potential for the Gibbs state is easily calculated if we 
explicitly write the grand partition function ( z = exp [(3[i\): 



oo N M M M 

7V=0 11=1*2=1 «JV = 1 



(N N 
a=l a<f} 

and observe that, because of the limitation of occupation to just the prime numbers (Eqj2]), 
the coordinates i a (a = 1, 2, • • • , N) can take values only among the 7r (M) prime numbers 
that exist in the interval [1, M\. Also, because of the impenetrability of the particles (EqJTJ), 
the choice can be made of 7r (M)!/ [tt (M) — N]\ manners. We obtain 

7r(M) 

z N = {l + zY {M) (6) 



w— n V / 



iV=0 

and 



oh] = -^iog(i+^r (M) . (?) 
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The prime-counting function 7r (x) is the number of prime numbers less or equal than x. 
In his 1859 classic memoir to the Berlin Academy of Sciences (see for example the monograph 



by Edwards|3] for a translation), Riemann gave an analytical formula for n (x) in terms of 
the zeros and the pole of £ (s) := X/n=i n ~ s i the zeta function named after him. Restricting 
to values of x such that x > 2 n with n = 1, 2, .. , it can be written |3j: 



7T (#) = 7T (x) +7T (x) , (8) 

where the smooth part, to which contribute the pole at s — 1 and the integer zeros of ( (s), 



is given by 



oo 

7f fx) 



n 

71=1 



(9) 



Li (zV») -^Li (x~ 2fe /") 
fc=i 

with Li(x) the logarithmic integral function and \i (n) the Mobius function. The complex 



zeros, on the other hand, contribute to the oscillatory part 



5? (x) = -2Re J] ^ ^ Li . (10) 

n=l a=l 

For x large enough the non-oscillatory and oscillatory parts can be written, respectively, 



as 



and 



7f(x)wli(x) (11) 



7T X 



-- — 22 ytw\ [°" q cos ^ q log x ) + ta sin ( ta log ^ ■ ( 12 ) 

iogx^ =i ^<t q + t a j 

For simplicity, from now on we will consider this last situation; so M will be assumed large 
enough. 

In Eqs. fllOp and (|T2|) the real numbers a a and t a are the real and imaginary part, re- 
spectively, of the complex number p a = o a + it a that verifies ( (p a ) = 0. The properties 
of the complex or non trivial zeros of the Riemann 's zeta function has been extensively 
studied jsj]. For example it has been demonstrated that there are infinitely many complex 
zeros and that all of them lie inside the region < 9? (s) < 1. Also is well known that if 
f (s) : = | ( s - 1) 7r- s / 2 r (§) C(s) then the non trivial zeros of ((s) are precisely the zeros of 
£ (s) and since £ (s) = £ (1 — s), we have that the complex zeros of £(s) are symmetric with 
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respect to the so called critical line 3? (s) = 1/2. The Riemann Hypothesis is the statement 
that all the complex zeros have their real part exactly on the critical line: a a = 1/2 Vet. Al- 
though billions of complex zeros have been numerically calculated [7] confirming all of them 
this conjecture, the full demonstration of the general validity of the Riemann Hypothesis 
remains still open. 

We wish to analyze the behavior of the equilibrium grand partition potential (EqJTJ), with 
7r (M) given by Eqs. (18UI2]) . as a function of the quantities a a G (0, 1) (a = 1, 2, • ■ ■ ). These 
numbers must be thought as possible values for the real part of the zeta function zeros and 
our goal would be to identify the true ones. 

We focus into a generic zero p 7 and observe that, since Riemann actually didn't know the 
reliable values of the complex zeros, Eqs. ffTUj) and f[T2"j) must be taken as general enough as 
to contemplate the possibility that, associated to £ 7 = 3 (p 7 ), both cases, er 7 ^ 1 — a 7 and 
<7 7 = 1 — <7 y , can be considered for (p 7 ). 

Explicitly we rewrite Eq.flj[J) for 0,M and z fixed 

09, M, z; <r 7 ) = ~ log (1 + zf M ^ , (13) 

with 



7T (M; o- 7 ) = ?f (M) + tt^ 7 (M) + ^ 7 (M; <r 7 ) . (14) 

Here the term 7r^ 7 (M) includes all the complex zeros except that labelled 7 and we assume 
that its contribution to the prime-counting function is the true one. We write the remaining 
term: 



2 / M a -i 



logM 



vr 7 (M; (T 7 ) w <( n M(l v ?L [(1 - <r 7 ) cos (t 7 log M) + t 7 sin (t 7 log Af )] ) for a 7 ^ 1 - a. 



[cr 7 cos (t 7 log M) + t 7 sin (t 7 log M )] + 



(l-(T 7 )^+t 



2 Af 1 / 2 ri 



log Af 



Tj-p [| cos (t 7 log M) + t 7 sin (t 7 log M)] for cx 7 = 1 — cr 7 , 



(15) 

where the fact that if ( (<r 7 + zt 7 ) = then also is ( (1 — ct 7 + ii 7 ) = 0, has been considered. 
Special cases can make the analysis even simpler. For M = Mi = expf^ {An — 1) / |t 7 |] 
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and M = M 2 = exp[| (4n — 3) / |£ 7 |] (n e N large enough) and taking into account that 
\t n \ > 14.37 for all the non trivial zeros of the zeta function, we have 



(m% + M^~ a A for <7 7 ^ 1 - <7 7 



7r 7 (M 1)2 ;(7 7 ) « ^ 2 1/2 ^ ^ 

fc |t 7 |logM 1|2 M l,2 fOT ff 7 = 1 " ^ 

where the upper sign corresponds to M\. Note the avoidable discontinuity in the function 
7r 7 (M; cr 7 ) when cr 7 = 1 — a 7 . 

Using Eq. ffTBl together with Eq. ffTol) -or the particular cases given by Eq. tfTB]) - is easy to 
see that the grand potential has a unique extremum (minimum or maximum depending on 
M) at the interval < er 7 < 1. From a heuristic point of view it is reasonable to expect 
that this extremum (maximum = unstable equilibrium; minimum = stable equilibrium) be 
reached when <r 7 takes the value (c" 7 ) <rue that gives, through Eq. (|T4|) . the correct number of 
prime numbers lying inside the interval [1,M]. The same conclusion can be achieved in a 
more formal way if we take into account that fixing a given value for <r 7 can be thought as 
imposing a constraint that limits (through 7r (M; <t 7 )) the number of the accessible states to 
the system. If this constraint is removed, so a 7 is leaved to freely vary in the interval (0, 1), 
then the equilibrium probability -P(er 7 ) that the system be in states with the parameter 



taking values in the interval between a 7 and <r 7 + 5a 7 behaves as|8] 

f exp[-/3ft09,M,z;<7 7 )] if a 7 e (0, 1) 
PK)0C \ ifa 7 ^(0,l). 

Moreover, cx 7 can be treated as a stochastic variable a 1 = a 7 (t) describing, for example, 
the position of an hypothetical Brownian particle of mass m moving in a force field with 
potential function given by 

/ Sl(J3,M,z;<r y ) if a 7 G (0, 1) 
^^^l oo ifa 7 ^(0,l). (18) 

In this picture, -P(cx 7 ) is the marginal of the joint probability distribution 

f [ 1.2 

P(<T 7 , <J 7 ) = if exp < —f3 U (<j 7 ) + -ma^ 



(19) 



where <r 7 = —2- and if is the normalization constant. By direct substitution one can verify 



that this probability distribution is the stationary solution of the Fokker-Planck equation 



9] 
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d_ 

at' 



P ( cr T , cr 7 , 1 1 = £ ( cr 7 , <r 7 , 9-, o> • P ct 7 , cr 7 , t 
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1 9 



= < -cr. 



7 ^7 m ^ 



da~y 



+ Ao% 



+ 



f 3 ™ 2 da' 



P(a 7 ,a 7 ,t), (20) 



where A is a friction coefficient. 

Associated to the Fokker-Planck equation (1201 with the initial condition P(cr 7 , cr 7 ,t ) 
5 (er 7 — <r 70 ) <5(cr 7 — <7 70 ) we have[l0j the Langevin stochastic differential equation 



m- 



<ia 7 



dC/ (a 7 
da-y 



^y+\ljf(t) 



(21) 



where / (i) is an additive Gaussian stochastic process <5-correlationated and of zero mean. 
Here — Ao~ 7 represents the friction force that affects the Brownian particle. Given the arbi- 
trariness (to our purpose) of the friction parameter, we simplify the analysis by choosing as 
such A = so that Eq. (|2ip yields the simpler Newton deterministic equation 



m- 



d(Ty 

~dt 



dU (<r 7 ) 
da^ 



(22) 



If we assume that initially the particle is placed at the position a 70 = (c 7 ) ea . tr that extremises 
Q (/3, M, z; <t 7 ) with velocity cr 70 = 0, then the equation of motion has as the unique obvious 
solution the isolated fixed point <r 7 (t) = (o" 7 ) . Vt. Taking into account that for any other 
pair of initial conditions the solution is time-dependent and that in the mathematical world 
the zeros of ( (s) do not change with time but remain constant forever, we deduce that it 



must be (a, 



7 7 true 



(cO . . One then infers that (cr~). should verify 

\ 7 / extr \ i 'true J 



dn M, z; a 7 



1 , ,„ N 97r 7 (M; cr 7 ) 



0. 



(23) 



= (°'7)t 



T 7 = (cr 7 ) true ^ 7 

The symmetry with respect to <r 7 and I — o~ 7 of expressions (|15|) or ffTBT) implies that this 
equation has the form [/ (<r 7 ) — / (1 — cr 7 )] , ^ =0 (with / (<r 7 ) ^ 0; / (1 — <7 7 ) 7^ 0) 
which is clearly verified by (o r 7 ) irue = 1 — (°"7) true — 1/2- Because the zero p 7 that we have 
considered is a generic one, we would conclude that (<J a ) true = 1/2 Vet, say that, in fact, the 
real part of all non trivial zeros of Riemann's zeta function lie on the critical line. 
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